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Résumé. — Cet article est consacré à l'étude des mesures limites associées à la solution de 
l'équation de Helmholtz avec un terme source se concentrant en un point. Le potentiel est supposé 
C°° et l'opérateur non-captif. La solution de l'équation de Schrôdinger semi-classique s'écrit alors 
micro-localement comme somme finie de distributions lagrangiennes. 

Sous une hypothèse géométrique, qui généralise l'hypothèse du viriel, on en déduit que la 
mesure limite existe et qu'elle vérifie des propriétés standard. Enfin, on donne un exemple 
d'opérateur qui ne vérifie pas l'hypothèse géométrique et pour lequel la mesure limite n'est 
pas unique. Le cas de deux termes sources est aussi traité. 



1. Introduction. 

Dans ce papier, on étudie la limite hautes fréquences de l'équation de Helmholtz, en dimen- 
sion n > 1, avec un terme source se concentrant en x = 0. Par un changement de variables, 
ce problème se ramène à l'étude de la limite quand h — > de Uh, la solution de 

(-!ÇA + V(x)-E)u h = h- n / 2 sQ. 

Le membre de droite de cette équation est une fonction qui joue le rôle d'un profil se concen- 
trant en à l'échelle h. Sa norme L 2 est préservée. L'énergie, ou paramètre de régularisation, 
E est dans le demi-plan complexe supérieur et vérifie E = Eq + hE\ + o(/i), Eq > 0, quand 
h — > 0. De plus le comportement du potentiel V est prescrit à l'infini (dans cet article, V sera 
à longue portée). 

Plus précisément, on veut calculer les éventuelles mesures semi-classiques de Wigner as- 
sociées à Uh- Autrement dit, on cherche à savoir si, pour q G C y3 (M? n ), la quantité 

lim h(Op(q)u h ,u h ), 

existe et s'écrit (q,fi h -i) ou ^h,- 1 es t une mesure sur T*M n = M. 2n indépendante de h (Op(q) 
est l'opérateur pseudo-différentiel défini en (jl.3p ). La présence du h devant le produit scalaire 
sera obligatoire pour rendre la quantité convergente. 



Classification mathématique par sujets (2000). — 35Q40, 35JfO, 35S30, 8fQ20. 
Mots clefs. — Equation de Helmholtz, Mesure semi-classique, Analyse micro-locale. 
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En général cette mesure vérifie les propriétés suivantes. Elle est supportée dans la surface 
d'énergie Eq. Elle est nulle dans la zone entrante (on parle aussi de condition de radiation à 
l'infini). Enfin cette mesure limite vérifie l'équation de Liouville : 

(£3* - VV{x)d i +21111^)^-1 = {2TT) 1 - n \S{O\ 2 5 x=0 5 l _ e+vm=Eo . 

On verra que le dernier point n'est pas toujours satisfait. 

A notre connaissance, le premier travail dans cette direction est du à Benamou, Castella, 
Katsaounis et Perthame pp. Sous l'hypothèse que V = 0, ils obtiennent la convergence vers 
la mesure précédemment décrite. Ce résultat a été étendu au cas des potentiels variables par 
Wang [25j . Lorsque la source est portée par une variété, des résultats équivalents ont été 
obtenus par Castella, Perthame et Runborg |6j dans le cas V = et généralisés par Wang et 
Zhang [26j aux cas à coefficients variables. Enfin, quand le potentiel est constant de chaque 
coté d'un hyperplan sur lequel il présente une discontinuité, la propagation de telles mesures 
semi-classiques a été traitée par Fouassier [10j . Tous ces travaux utilisent des majorations de 
la résolvante entre espaces de Besov (estimations de Morrey-Campanato) . Ces estimations 
sont parfois obtenues à l'aide de la théorie de Mourre. C'est pourquoi la régularité du potentiel 
peut être abaissée jusqu'à C 2 , dans certains cas. Néanmoins, tous supposent que le potentiel 
est nul ou satisfait l'hypothèse du viriel. Ce présent travail a pour but de décrire ce qui se 
passe quand cette hypothèse est retirée, mais que l'opérateur reste non captif. 

Mentionnons qu'il existe également des articles qui prouvent spécifiquement des estimations 
de la résolvante entre espaces de Besov (voir, par exemple, Perthame et Vega |20] pour des 
potentiels ne s'annulant pas à l'infini, Wang et Zhang [26J pour des termes sources portés par 
des surfaces, Castella et Jecko \4\ dans le cas C 2 non-captif, Wang |25] et Castella, Jecko et 
Knauf [5j pour des singularités coulombiennes) . Enfin, Castella [3] a étudié la limite faible 
hautes fréquences de l'équation de Helmholtz dans le cas non captif. 

Dans cet article, on démontre que la solution de l'équation de Schrôdinger stationnaire peut 
s'écrire micro-localement près de tout point po £ ^ 2n comme somme finie de distributions 
lagrangiennes. Toutefois, la variété x = pose des difficultés, et on peut juste majorer la solu- 
tion près de l'origine. Une fois cette construction effectuée, il est facile de déduire l'existence 
de la mesure limite et de la calculer. A cause du caractère bilinéaire de ces mesures, on doit 
faire une hypothèse géométrique sur les variétés lagrangiennes qui portent la solution. Cette 
hypothèse est plus générale que l'hypothèse du viriel. Elle, ou une autre du même type, est 
en fait obligatoire : dans la partie HJ en s'inspirant de l'exemple de Castella PU section 9], on 
construit un opérateur qui ne vérifie pas l'hypothèse géométrique et pour lequel la mesure 
limite n'est pas unique. 

Pour obtenir ces résultats, on utilise une méthode dépendante du temps et on écrit la 
résolvante comme l'intégrale du propagateur en temps positifs. Des constructions BKW 
développées par Maslov et Hôrmander servent alors à approcher le propagateur. Cette idée 
avait été utilisée par Castella dans [3], mais elle vient surtout, pour ce type de questions, de 
l'étude de l'amplitude de diffusion réalisée par Robert et Tamura |24j en 1989. Remarquons 
que, en un certain sens, l'amplitude de diffusion peut être vue comme l'équation de Helmholtz 
avec terme source à l'infini. Ces méthodes, très gourmandes en dérivées, nécessitent que le 
potentiel et la transformée de Fourier du terme source soient essentiellement C°° . Par contre, 
la preuve n'utilise ni la théorie de Mourre ni les estimations de la résolvante entre espaces de 
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Besov. Cependant le lemme [375], qui sert à contrôler la norme 1? (d'une partie) de la solution 
très près de x = utilise une décomposition dyadique. 

Enfin, avec cette méthode, on peut aussi traiter le cas de deux (ou plusieurs) points sources. 
Nous aurons besoin d'une hypothèse géométrique qui, comme précédemment, est essentiel- 
lement nécessaire pour l'unicité de la mesure. Fouassier [9] avait étudié cette question sous 
l'hypothèse V = 0. 

On se place dans le cadre de l'analyse semi-classique et h e]0, ho[ désigne le petit paramètre. 
Dans la suite, on utilisera les notations suivantes. Pour m(x, £) > une fonction de poids, 
S(m) désigne l'ensemble des fonctions q(x,Ç,h) G C°°(R^) telles que 

(1.1) d^q(x,^,h)=0(m), 

pour tous a, (3 £ N n . On dit que le symbole q(x, £, h) est classique et on note q G S c \(m) si il 
existe une suite de symboles qj(x,Ç) G S(m) tel que 

N 

(1.2) q(x,Ç,h)-J2<îj{ x ,0h j G S{h N+l m) :=h N+l S{m). 

j=0 

Sous cette hypothèse, qg est appelé le symbole principal. Pour q G S (m), l'opérateur pseudo- 
différentiel (semi-classique en quantification de Weyl) de symbole q est défini par 

(1.3) Op(q)f(x) = 1 1 e^-y^ h q(^,i,h)f{y)dydi, 

au sens des intégrales oscillantes. On note ^(m) = Op(5(m)), l'ensemble des opérateurs 
pseudo-différentiels de symbole dans S (m). Enfin, 

(!- 4 ) m = ^Je-^ h f(x)dx, 

est la transformée de Fourier de / G 5'(lR ra ), le dual de l'espace de Schwartz 5(M n ). 



2. Enoncé des résultats. 

On considère un opérateur de Schrôdinger semi-classique sur L 2 (M ra ), avec n > 1, 

(2.1) P = _^A + y(x), 

et on fait les hypothèses suivantes. 

(Hl) V est une fonction C°° sur M. n , et il existe p > tel que 

|0°V(z)| < (x)""-'"!, 

pour tout a G N n . 

Autrement dit le potentiel est régulier et à longue portée. Le symbole de P est noté 

p(x,0 = Ç + V(x), 
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et son champ hamiltonien 

(H2) La surface d'énergie Eq > est non captive. C'est à dire, 

)C(Eq) = {(x,£) G p~ 1 (E ); t — > exp(iH p )(x, £) reste borné en temps} 

(H3) De plus, V(0) < E . 
(H4) On suppose que 



mes n _i{£ G v / 2(^r r V(Ôj)S"~ 1 ; 3t > U x exp(iH p )(0, £) = 0} = 0. 

mes n _i est la mesure de Lebesgue sur y/2(Eo — V(0))§ n_1 et U x (x,^) = x est la projec- 
tion spatiale. 

Notons que l'ensemble qui apparait dans |(H4)| est mesurable car fermé. L'interprétation 
géométrique de cette hypothèse est la suivante (voir la partie [6] pour les preuves). La so- 
lution de l'équation de Helmholtz sera micro-localisée sur l'ensemble A défini par 

(2.2) A := {exp(ffl p )(0,0; t > et Ç 2 = 2(E - V(0))} C p" 1 ^)- 

Cet ensemble peut se "recouper" (voir les figures HJ [3] (en dimension n = 1) et [5] (en 
dimension supérieure)). Mais localement près d'un point pç, = (xo,£o) e P^H-^o)) il est de la 
forme 

K K 

(2.3) A = (J A fc , si x ± et A = X U |J A fc , si x = 0, 

k=l k=l 

où les ensembles A et A& sont définis ci-dessous. La variété 

(2.4) À = {exp(tH p )((U); < t < e et Ç 2 = 2(£ - K(0))}, 

avec eo > assez petit, est une variété (locale) lagrangienne C°° se projetant bien suivant les 
variables d'espace, pour x ^ 0. Elle est représentée sur la figure CD Les A& sont des variétés 
(locales) lagrangiennes C°°. Plus précisément, 

(2.5) < t, Ç 2 = 2(£ - V(0)) et exp(tH p )(0, Ç) = p }, 
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Figure 2. Un exemple de potentiel V qui, pour < Eq < max V, ne vérifie pas 
l'hypothèse du viriel mais pour lequel [(H4)| est vrai. En effet, l'ensemble qui apparait 
dans |(H4)| est un singleton {£o}- 



est un ensemble fini {(ti, . . . , {tK-, (,k)} et 

(2.6) Afe = {exp(tH p )(0, £); t au voisinage de tk et £ au voisinage de Ç^}. 

Ainsi, il se peut que Aj = A& avec j ^ k. On notera parfois Ao = A et, d'après la section [6j 
l'hypothèse |(H4)| est équivalente à 

(H4)' Les intersections des variétés A et A^ sont de mesure nulle : 

mes n (Aj n A fc ) = 0, 
pour < j k < K. La mesure est prise sur la variété Aj (ou A*.). 

D'après le principe d'absorption limite, on sait que, pour tout E e]0, +oo[, 

(2.7) (P - {E + iO))- 1 := lim(P - {E + itf)) -1 , 

existe en tant qu'opérateur de L 2 ((x) 2a dx) dans L 2 ((x)~ 2a dx), pour tout a > 1/2. De plus, 
Robert et Tamura |22| ont démontré que, sous les hypothèses piTJI et [(H2)| 

(2.8) sup ||(x)- a (P- (P + ^))- 1 (x)- Q || <h-\ 

0<S<1 

pour E g]0, +oo[ au voisinage de Eq. Leur preuve utilise, entre autres choses, les constructions 
d'Isozaki et Kitada. Mais il existe d'autres démonstrations, notamment celle de C. Gérard 
et Martinez |12| qui est basée sur la méthode de Mourre [19] . D'autre part, Burq [2] a 
prouvé (|2.8p par l'absurde grâce aux mesures de défauts. Cette estimation résulte aussi du 
récent papier de C. Gérard qui en donne une preuve directe utilisant uniquement de 
l'analyse fonctionnelle. Dans la suite et dans ce genre d'estimation, a désignera toujours un 
réel strictement supérieur à 1/2. 

On prend une donnée Sh qui se concentre en de la manière suivante : 
(H5) On considère la suite de fonctions 

S h (x) = h- n ' 2 s(l), 

où S est dans {x)- N L 2 (R n ) pour tout N > 0. Donc, S G C°°(R n ). 
et on fait l'hypothèse suivante sur l'énergie : 
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(H6) Soit E = E(h) £ C une famille d'énergies vérifiant 

E = E + hE x + o(h), 

et ImE > 0. 

Ce papier est consacré à l'étude de la solution de Helmholtz 
(2.9) u h (x) = (P-(E + iO))- 1 S h . 

Si ImE > 0, on a u h = (P - E)~ l S h . 



Théorème 2.1. — Sous les hypothèses (H1)-(H6) il existe une mesure de Radon positive 



fj, h -i sur M. 2n telle que, pour tout q G C£°(R* n ), 

(2.10) h(Op(q)u h ,u h ) — > (q^h- 1 ), 

quand h —>■ 0. La mesure u h -i ne dépend pas de h et l'indice est juste une notation. 
La mesure p h -i possède les propriétés suivantes qui la caractérisent uniquement : 

i) La mesure u h -i est supportée dans la surface d'énergie Eç, : 

(2.11) suppfi h -i Cp^iEo). 
H) La mesure /x^-i est nulle dans la zone entrante : 

(2.12) n h ~i = dans {(x,Ç); \x\ grand et cos(x,Ç) < —1/2}. 
iii) La mesure fi h -\ vérifie l'équation d'évolution : 

(H p + 2 Im Ei) n h -i ={id x - VV{x)d i + 2 Im Ei)fi h -i 



(2.13) =(2vr) 1 -"|5(Ç)|^ :E= o<5i ç2+y(0) 



--E - 



Le membre de gauche de l'équation (|2.10p est à prendre au sens de la dualité L 2 ((x) 2a dx), 
L 2 ((x)~ 2a dx). On dira qu'une fonction est un C Q (1) si sa norme L 2 ((x)~ 2a dx) est un 0(1). La 
mesure [i h -i est indépendante de h. L'indice sert juste à indiquer que le produit scalaire 
dans le membre de gauche de ()2.10j) a été multiplié par h. 

Remarque 2.2. — En fait, on calcule micro-localement Uh près de tout point po = (xq, £q) e 
R 2n avec xq ^ (pour xq = 0, on a juste les majorations données par le lemme 13.51 et la 
discussion après le lemme I3J1 (voir aussi la proposition \2.3\) ). Si p(po) = Eq et xq / 0, aiors 
Uh est somme de K distributions lagrangiennes d'ordre h~ï , à symbole classique, portées par 
les variétés lagrangiennes (voir (|3.47[) et (|3.60p )■ D'autre part, Uh est nul micro-localement 
près des points po £ p~ 1 (Eo) avec xq ^ 0. Remarquons que ces résultats ne nécessitent pas 
l'hypothèse \(m)\ 

Comme corolaire, ces constructions donnent un développement asymptotique (au sens C°°) 
explicite de Uh si les lagrangiennes se projettent bien spatialement. Enfin, la preuve du 
lemme 1331 doit permettre d'étendre la construction de la distribution lagrangienne portée 
par A jusqu'à des voisinages de taille h 1-6 de faisant apparaître une phase singulière en 
comme dans le cas de —h 2 A/2. 
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L'hypothèse (H4) , ou une hypothèse du même type, est obligatoire pour obtenir les conclu- 



sions du théorème 12.11 L'exemple donné dans la partie [4] montre que, si cette hypothèse est 
retirée, la mesure limite n'est plus unique. 

Dans la partie consacrée à l'étude des mesures semi-classiques de l'équation de Schrôdinger 
stationnaire de son papier |25j . Wang démontre le théorème 12. Il en supposant que le potentiel 
V G C 2 (R n ) satisfait 

(2.14) u (E - V(x)) - x ■ W(x) > c , 

avec c > et u G]0,2] (équation (4.23) de [25]). Comme p = £ 2 + V(x) dans [25J, si 
(x(t),£(t)) G p~ 1 {Eq) est une courbe hamiltonienne, il vient 

d t (x(t) ■ f(t)) =2Ç(t) 2 -x-VV(x) 

>vot(t) 2 ~ x ■ W(x) 

(2.15) =v Q (E Q -V(x))-x - VV(a?) > c . 

En particulier, x ■ £ est strictement croissant en temps et l'ensemble qui apparait dans l'hy- 
pothèse [(H4j| est vide. 

L'hypothèse | (H4)] est proche de l'hypothèse de transversalité énoncée par Castella dans [3]. 
Cet article est consacré à la limite faible de l'équation de Helmholtz ; autrement dit, à la limite 
de {uh,h~ n / 2 ip(f^)), pour toute fonction tp G C^°(M n ). L'auteur suppose que la dimension n 
est supérieure à 2, que l'ensemble 

(2.16) M := {(t, n); Ç 2 = 2(E - V(0)) et exp(tH p )(0, = (0, rj)}, 

est une sous-variété de M 2n+1 de codimension strictement supérieure à n + 2 et que l'espace 
tangent de A4 est donné par la linéarisation des équations qui le définissent (hypothèse (H) 
page 245 de |3j). Mais alors, l'intersection de deux variétés Aj et A^, pour < j / k < K, est 



une sous- variété de Aj (ou A&) de codimension supérieure à 1 et (H4)' est vérifiée. Néanmoins, 
le problème traité par Castella est différent de celui abordé ici et les comparaisons ne sont 
pas évidentes. 

Le théorème 12.11 décrit la contribution principale de la solution de l'équation de Helmholtz. 
Mais il reste une partie sous-principale qui n'est pas principalement portée par la surface 
d'énergie Eq. Cette contribution est décrite par la proposition suivante. 



Proposition 2.3. — Sous les hypothèses (Hl) - (H6) , il existe une mesure de Radon positive 
/il sur R 2n telle que, pour tout q G (^(M 2 ™), 

(2.17) (Op(?(p - E ) 2 )u h , u h ) — > {g, 

quand h — > 0. De plus, on a 



(2-18) (g, Ml )= J q(0,0\S(H)\ 2 dti. 

Le fait que Uh porte une contribution particulière en x = est spécifique à l'étude de la 
résolvante. Dans le cas de la mesure spectrale, 

(2.19) V(E) =±-(( P -(E + îO))- 1 — (P — (E — iO))- 1 ) , 
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E G]0, +oo[ proche de Eq, ce terme n'existe pas et V(E) est un opérateur intégral de Fourier 
régulier. Ce résultat a été démontré par C. Gérard et Martinez dans [13] . 

De plus, il semble aussi possible de définir une mesure [i h -\/i par 
(2.20) h l l 2 (Op{q(p-E))u h ,u h ) — (g,/x fc -i/ a >, 

pour tout g G C^°(lR 2n ). Sous une hypothèse géométrique, cette mesure semble portée par 
p (Eq) n {x = 0}. En fait, cette mesure est un cas particulier de limite faible de l'équation 
de Helmholtz ; elle code le retour en de la variété A lorsque l'hypothèse de transversalité 
(H) de Castella [3] n'est pas vérifiée. 

Avec cette approche, on peut aussi traiter le cas de deux (ou plusieurs) termes sources 
concentrés en x\ ^ x 2 , deux points de R n . Pour S 1 , S 2 deux fonctions satisfaisant l'hypothèse 
|(H5)| on pose 




En plus de l'hypothèse |(H4)| aux points x\ et X2, on rajoute : 

(H7) On suppose que 

mes n _i{Ç G y/2(E - V(x 1 ))S n - 1 ; 3* > U x exp(tH p )(xi, £) = x 2 ) = 0. 

mes n _i est la mesure de Lebesgue sur y/2(Eo — V(xi))S n ~ 1 . 

Autrement dit on suppose que les trajectoires dansp -1 (.Eo) qui vont de x\ à x 2 sont de mesure 
nulle. L'hypothèse |(H7)| ne change pas si x\ et x 2 sont intervertis. 

Pour j = 1,2, soient A^, k = 0, . . . , Kj les variétés lagrangiennes définies en (|2.4p et (|2.6p 
en po avec point de base xj. L'hypothèse |(H7)] est alors équivalente à 

(H7)' Les intersections des variétés A], et A 2 ^ sont de mesure nulle : 

mes n (A^ nA 2 fe ) = 0, 
pour tout < kj < Kj, j = 1, 2. La mesure est prise sur la variété A^ (ou A 2 ^). 

Dans ce cas, la solution de Helmholtz est alors la fonction Uh solution de 

(2.21) u h (x) = (P-(E + ï 0))- 1 {S 1 h + S 2 ). 

On a alors le théorème suivant qui étend le résultat de Fouassier [5] aux cas à coefficients 
variables. 

Théorème 2-4- — On suppose |(H1)| |(H2)| |(H6)| et |(H7)| On suppose aussi que les hy- 
pothèses p3)| [(H4)1 et [(H5)1 sont vérifiées en x\ et en x 2 (autrement dit, en remplaçant le 
point de base par x\ et par x 2 ). Alors, pour tout q G C^°(M. 2n ), 

(2.22) h(Op(q)u h , u h ) — > (q, nl-i + /^-i) quand h -> 0. 

La mesure H J h -i est la mesure fournie par le théorème \2.1\ avec point de base Xj. 
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La preuve suit celle du théorème 12. Il et on omet les détails. L'hypothèse | (H7)] implique que les 
interactions entre les deux points sources sont négligeables. Comme pour le théorème 12.11 si 
cette hypothèse est retirée, on peut construire un exemple d'opérateur pour lequel la mesure 
n'est pas unique. Remarquons que ces hypothèses impliquent n > 2. 

Le reste du papier est consacré à la preuve du théorème 12.11 Dans la section m on donne 
un exemple où la mesure n'est pas unique. La proposition 12.31 est démontrée dans la section 
Enfin, l'équivalence entre [(H4)] et [(H4) ' est prouvée dans la section El 



3. Preuve du théorème 12.11 

Dans les sections l3.1H3.6| on pourra supposer que E = Eq + hF%, où F\ G {z G C; Imz > 
0} ne dépend pas de h. En effet, toutes les constructions seront régulières par rapport au 
paramètre F± et toutes les estimations seront uniformes en F\ , quand F\ varie dans un compact 
de {z G C; Imz > 0}. Et pour obtenir le cas général |(H6)[ il suffira de prendre F\ = 
h~ 1 {E - Eq) =Ei + o(l), où E est donné par |(H6)[ 

3.1. Localisation près de la surface d'énergie Eq. 

On note Pq = —h 2 A/2 le laplacien libre de symbole po(x,£) = £ 2 /2. 

Lemme 3.1. — Soit f,g G Cq°(R) tel que f = 1 près de E Q - V(0) et g = 1 près de E . 
Alors 

(3.1) u h (x) =(P ~{E + iO))- 1 f(P )S h + O a {l) 

(3.2) =g(P){P ~ (E + iO)r 1 f(P )S h + Û (1), 
pour tout a > 1/2. 

Démonstration. — Soit 1e -< / G Cq°(R) tel que f(po) = 1 près du support de f(p). Par le 
calcul fonctionnel et comme S h = 0(1), on a 

u h (x) ={P - (E + iQ)y l j{P)S h + (P- E)-\l - J(P))S h 

(3.3) =(P-(E + i0))~ l f(P)S h + O(l). 

Ici, on a utilisé le fait que f(P) G ^(1) et donc f(P)Sf l G L 2 ((x)^ 0O dx). D'après le calcul 
fonctionnel des opérateurs pseudo-différentiels ([7] par exemple), on a 

(3.4) f(P) G f(P)f(P ) + ^(er 00 ). 

En utilisant ([23]), {33J), dUH]) et encore <Ç£ty . il vient 

u h (x) ={P - (E + iO))- l f(P)S h + 0(1) 

= (P — (E + iO))~ 1 f(P)f(P )S h + 0(1) + O a (h°°) 

(3.5) =(P - (E + iO))- 1 f(P )S h + 0(1) + O a (h°°), 

ce qui prouve (|3.ip . L'identité (|3.2p découle de (|3.3p . □ 
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Remarquons que, comme (P — E) 1 (1 — f(P)) G 2 ), la preuve montre que 

(3.6) u h (x) = (P-(E + iO)y 1 f(P )S h + R, 

où R = O a (l) vérifie (1 - (p(x))R = O a {h°°) pour toute fonction ip G Cg°(R n ) telle que 
<p = 1 près de 0. En fait, comme (P — (E + îO)) -1 = 0(h~ v ) en tant qu'opérateur de 
H s ({x) 2a dx) dans H s {{x)~ 2a dx), pour tout s > 0, on peut montrer que d^R = O a {h~^) et 
(1 - <p(x))d£R = O a {h°°) pour tout (3 G N n . 

3.2. Contrôle des temps grands. 

Comme dans Robert-Tamura [24] et Castella [3], on voit la résolvante comme l'intégrale 
en temps du propagateur. En effet, 

(3.7) (P - (E + M))' 1 = 1 [ T é-^ p -^ h dt + (P-(E + iO))- 1 e-* r ( p - E V h , 

h Jo 

pour tout T G M. On a alors le lemme suivant. 
Lemme 3.2. — Pour tout q G C£°(M 2n ), 

(3.8) Op(ç)(P — (-E + iO))" 1 /^/, = Op(q)~ £ e'^'^ f(P )S h dt + _/3 (/i°°), 
pour tout T assez grand et P > 0. 

Démonstration. — La preuve est entièrement donnée au bas de la page 173 de |24| . et on 
rappelle juste les deux étapes. D'abord, on énonce une version du théorème d'Egorov. 

Proposition 3.3 (Proposition 3.1 de [24]). — Soit uj(x,£) G Cfi°(R 2n ). Supposons que 
u>t G 5(1) soit nul près de 

{(x,0'i ( x ,0 = exp(iH p )(2/,77) avec (y,n) G suppw}. 

Alors 

\\(xf Op{uj t )e- itp/h Op(u)(xf\\ = 0(h°°), 

pour tout > 0. Cette estimation reste uniformément vraie pour t dans un compact de M à 
condition que u>t soit uniformément dans 5(1). 

On dit qu'un symbole uj G C co (M 2n ) est de classe Aq si 

\^u(x,0\<(x)-^(0- L , 

pour tout L > 0. Les ensembles T_ (resp. T + ) définis par 

T±(R, d, a) = {(x, £); |x| > R, cT 1 < \Ç\ < d et ± cos(x, £) > cr}, 

sont appelés les zones entrantes (resp. sortantes). Dans [24], la proposition 13.31 est énoncée 
pour uj t G -Ao, mais elle reste vraie pour oj t G 5(1). 

Comme les surfaces d'énergie proche de Eq sont non-captives, cette proposition implique 
que le micro- support de e- iT ( p - E V h f(P )S h (x) est micro-localisé dans une zone sortante 
modulo 0~P(h°°), pour tout f3 > 0. Le lemme 1331 est alors une conséquence directe du lemme 
suivant qui repose sur les constructions d'Isozaki et Kitada. 
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Proposition 3-4 (Lemme 2.3 de |24j). — Soit lu, uj± G Aq tel que le support de u>± (resp. 
lu) soit dans T±(R,d,a±) (resp. {\x\ < 9-R/10}), avec R assez grand et cr_ < o + . Alors 

\\(x)P Op(u T ){P - (E±iO))- 1 Op(u±){x) f3 \\ = 0(h°°), 

||(x) /3 Op(cu)(P- {E±iO))- 1 Op(io ± )(x)> 3 \\ = 0{h°°), 

pour tout (3 > 0. 

Remarquons que Robert et Tamura suppose que l'énergie E est réelle et que le potentiel est à 
courte portée, p > 1 (pour pouvoir développer la théorie de la diffusion), mais les propositions 
précédentes sont encore valables pour des potentiels à longue portée (voir l'appendice de |23j 
où la construction des propagateurs approchés entrant et sortant est faites pour p > (voir 
aussi les travaux d'Isozaki et Kitada |16j et |17j )) et pour des énergies vérifiant |(H6)| (en fait, 
pour E G C proche de Eq et ± ImE > 0). □ 

3.3. Contrôle des temps petits. 

En utilisant les lemmes [331 et 1331 et les propriétés du calcul pseudo-différentiel, on obtient, 
pour q G C^°(E. 2n ), 

(3.9) V (q)u h = Op(q)ù h + 0~ p (l), 
où ûh est défini par 

(3.10) ù h := ~ £ e~^ p - E V h f{P )S h dt. 
Dans cette section, on étudie 

(3.11) B 1 (x) :=^j\(t)e-^ p - E ^ h f(P )S h dt, 
avec x G C§°(M). 

Remarquons que, pour toute fonction g G C^°(IR n ) telle que g = 1 près de 0, 
f(Po)S h =g(x)f(P )S h + 0(h co ) 

(3.12) =(2vr)- n /i-5 j g{x)e ix < /h f(f/2)S(0d^ + O{h c °). 

De plus, le support de g (resp. /) est aussi proche de (resp. Eq — V(0)) que l'on veut. Donc 



ih- 1 - 



2 



T 



(3.13) Bl {x) = ^-^- x(t)e-^ p - E ^ h é x <l h g{x)f(e/2)S{i)dtdi + 0{h^). 

La norme 1? de B\ très près de x = est contrôlée par le lemme suivant. 

Lemme 3.5. — Soit x G Co°(^); X = 1 P r ^s de et de support suffisamment proche de 0. 
Alors pour tout e > 0, 

(3.14) l [x{<e B 1 = 0(<Sëh-$), 
en norme L 2 . 
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Démonstration. — D'après la méthode BKW (voir [8] ou [7] par exemple), on sait que, pour 
tout t assez petit, 

(3.15) e -it(P-E)/h ( g ( y yy</ h ) = a(t, x, Ç, h)é^ x ^l h + 0{h°°), 

La phase (p et le symbole a sont C°° par rapport à t, x, £. La phase vérifie l'équation eikonale 

j dt<p + p(x, d x tp) - E = Q 
\ P(0,a;,£) = x-Ç 

L'amplitude a(t,x,Ç) est à support compact en x et développable par rapport à h 



(3.16) 



a(t,x,Ç,h) ~ ^Vo,(t,:z,£). 

3=0 

De plus, les coefficients ctj vérifient les équations de transport habituelles 

(d t + V x <p ■ V x + A x ip/2 - Fi)ao =0 

(3.17) 1 

(5* + Vs;¥? • V x . + A x ip/2 - Fi)aj =-Aoj_i pour j > 1, 

avec données en t = 

(3.18) °o(0, 

aj(0, x) =0 pour j > 1. 

En particulier, comme annoncé au début de la partie O les a,- sont holomorphes par rapport 
à F x . 

On est alors amené à calculer 

(3.19) J := [T x(t)e^' fl/h a(t, x, £, fc)/(É 2 /2)S(0 d£ dt, 



J Jo 

— l— n 

puisque B\ = %h ^Ç J + 0(h°°). Pour ô > h petit, on note 

(3.20) j s (x) := l 5<lxl<2S j£ x(t)e l ^ x ^ h a(t,x,Ç,h)f(e/2)S(0dÇdt, 

et on essaye d'estimer J$. 



Bien sûr, pour x £ Cg° 



T 



(3.21) Jf(x):=l w<5 x(t/%W^ (M ' â/ "a(^,Ç^)/K 2 /2)S(Ç)^rft, 
vérifie 

(3.22) \\Jf°\\^=0{S 1+ ^). 

Remarquons que, d'après (|3.16|) . on a 

(3.23) <p(t, x,Ç)=x-Ç- t(f/2 + V(x) - E ) + t 2 S(l). 

De plus, grâce au lemme [37T| Ç 2 reste dans un compact disjoint de sur le support de /(£ 2 /2). 
En particulier, pour C\5 < t < 1/C% avec C\ assez grand (indépendant de 5), on a 

(3.24) - £ • dçtp = -x • Ç + tf + t 2 S(l) > 5, 
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sur le support de ls<\x\<2â a {t> x >£) ^)/(Ç 2 / 2 )- Comme les dérivées par rapport à Ç de la fonc- 
tion son * uniformément bornées, des intégrations par partie en £ donnent, pour 
X G Cq°(] — oo, Ci + 1]) avec x = 1 près de ] — oo, C%], 



J+(x) :=l s<M<26 X(t/ô))x(t)e^ x ^ /h a(t, x, Ç, h)f(e/2)S(Ç) d£ dt 

(3.25) =0((h/ô) M ), 
pour tout M G N. Ainsi 

(3.26) \\4\\ L2 =0(ô^(h/ô) M ), 
pour tout M G N. 

D'un autre coté, pour t < 6/C2 avec C2 assez grand, on a 



• dç<p x 2 — tx ■ £ + t 2 x5(l) 



( 3 - 27 ) — /F = 39 — > 1. 



.T 2 



sur le support de l<5<|:r|< 2( 5a(t, Ci M/(£ 2 / 2 )- Comme les dérivées par rapport à £ de la fonc- 
tion x ï_ tx .ç+ t ï x s(i) sori t uniformément bornées, des intégrations par partie en £ donnent, pour 
X G Cq°(] - 00, C 2 ]) avec % = 1 près de ] - 00, C 2 /2], 



(3.28) =((W M ), 
pour tout M G N. Ainsi 

(3.29) \\Ji\\ L2 = 0[ô^(h/ô) M ), 
pour tout M G N. 

Il reste à étudier 

(3.30) Jl(x) := 1 5<W<25 j£ x(t/S) X (t)e l ^ x ^ h a(t, x, £, fr)/(£ 2 /2)£(£) d£ di, 

où x = (l ~~ x)x à son support compact disjoint de 0. On fait alors le changement de variables 

t = ÔT, x = ÔX, £ = H. 

Le symbole ls<\x\<2SX(t/S)x(t)a(t, x, £, /i)/(£ 2 /2),S(£) devient l 1<: \ x \ <2 â(T,X,E) où le sym- 
bole ô G S(l) est à support compact dans ]0, +oo[rx(IR n \ {0})x x (K™ \ {0})=- On pose 

(3.31) $(T, X, E, ô) = ô^ipiôT, ÔX, E). 

Cette fonction est C°° par rapport aux variables T, X, E et 5. D'après (|3.23j) . elle s'écrit 

(3.32) $(T, X, H, S) = X ■ E - T(E 2 /2 + V(0) - E Q ) + 5R(T, X, E, ô), 

où R G 5(1) sur le support de â. Toujours sur le support de a, les points critiques en (T,E) 
de <£> sont les points tels que 

/ - ~ 2 /2 - V(0) + E + ôd T R\ 
3.33 d T ~<S>=[ y =0. 

y ' ' \ X — TE + ôd^R i 
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Pour 5 = 0, il existe un unique point critique \T(X,0) = -j=J=B== ; S(X, 0) 

— * — I (car i > U sur supp(aj). h,n ce point, 

/ -*E 

(3.34) HessT = * = I _ c _ T 

est de déterminant — T"™ -1 - 2 C e qui est uniformément minoré sur supp(a). Par le théorème 
des fonctions implicite, pour 5 assez petit, l'équation dr,s^ = a donc une unique solution 
(T,È)(X,S) = (f,H)(X,0) + ÔR avec R G 5(1). En particulier, d'après <^32\i et (ET34"1) . 
HessT,s *£* est uniformément minoré au point critique pour (5^1. En écrivant 

(3.35) <9 T =$(T, S) = Hess T ,~ S(T, H)(T - f, E - S) + 0(|T - f | 2 , |3 - S| 2 ), 
on constate que 

|(r,H)-(f,s)| 



(3 - 36) \d Tr ^(T,~)\ ' 

est uniformément borné. En appliquant la méthode de la phase stationnaire (le paramètre 
semi-classique est ici h/S) k paramètre (théorème 7.7.5 de |14j ), il vient 

J${x) =ôjj l 1<|x|<2 â(T,X,H)e^^- 5 ) 5 /^Hdr 

(3.37) =0(ô(h/ô)^). 
Et donc 

(3.38) Ps^Wl* = 0(ô^ô(h/ô)^) = 0{5h s ^). 

En, utilisant une décomposition dyadique ô = 2~ n avec h 1 ^^ < 5 < 2e ainsi que (|3. 22p . 
(^26]) . dS^SD et (1338]) . il vient 

I|1n<^IU 2 <(W-mII + K + wJI) + ( E rai + h j °h + il J /n) 

5=2-" 
ft 1 -' l <5<2e 

<fc(l-M)(*+l) + ^(l^)I^M + £ 5 Ç (/i/J) M + + ^ (/l/(5) M 

5=2'" 
/i 1 -''<5<2e 

(/iè-Mf +D + h ^M-l-^ + h ,M-^ llQhl +VI+h ,M-^ llnhl) 

(3.39) ^h^y/i+h"), 

avec ^ > 0. Pour obtenir (|3.39j) . on fixe fi > assez petit (pour majorer le premier terme) 
puis M assez grand. □ 

3.4. Calcul près de 0. 



Dans cette partie, on reprend le calcul de 

(2tt)- 



(3.40) Bl {x) = 1 ——^ [[ X (t)e^ t ^ h a(t,x^,h)f(e/2)S(0dCdt + O(h o °), 
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pour ôi < \x\ < 82 avec < ô\ < 62 petits mais fixés. D'après (|3.29p . on a vu que 



• 7 — 1— - p pT 

ih a ' ' 



(;>K " 1 BlM) 1 1 \(t)e^'«/ ft a(t,x,e, h)f(e/2)S(0 < ,11 ! 0(1^ ). 







avec x = x(l ~~ x)> ou X £ Co°0^) et x = 1 près de 0. D'après (|3.30p ^ (|3.37p . 011 °l ue l' on 
peut utiliser la méthode de la phase stationnaire en t, Ç dans (|3.4ip (pour tout x, la phase a 
un unique point critique qui est non dégénéré). D'un autre coté, en un tel point, on a 



(3 42) f^(t,x,6 = 



Rappelons que, d'après la proposition IV-14 i) de [21] . 

(3.43) (x, d x (p) = exp(tH p )(ôçc/?, £), 
et d'après (40) de |21j . on a 

(3.44) Ç(a) = ^(a,x(a),0, 

pour toute courbe hamiltonienne (x(s), £(s)) avec £(0) = £. En utilisant (|3.16|) . (|3.43p et 
(|3.44p . on voit que (|3.42p est équivalente à 



(3.45) 



p(0,0 = E 
n x exp(tH p )(0,0 =x. 



En utilisant (|3.44j) et en notant 7(2;) = (x(s), £(s)) la courbe hamiltonienne exp(sH p )(0, 
donnée par (|3,45p . il vient 

tp(x) :=<p(t, x, = (p(0, x(0), + / d x (p(s, x(s),Ç)x(s) ds 

Jo 

(3.46) = / Ç(s)x(s)ds= I £dx. 

JO J-y(x) 

Et donc, la méthode de la phase stationnaire dans (|3.41|) donne 

(3.47) B^x) = b(x, h)h-\ê^l h + 0(h°°), 

où b(x,h) ~ J2j>obj(%)hi est un symbole classique à support compact avec 

i ï sgn w'! c 

(3.48) 6 = i(2vr) — £(t) Qo (t 5 x , 0/(^/2)5(0, 

| det ç| 2 

où (t,0 es t l'unique couple de points (avec t dans le support de x) vérifiant (|3.45p . A l'aide 
de (|3.34p . on trouve sgn<^^ = 1 — n et on remarque que 

1 n—1 

(3.49) |det^|3 =t~ |e|(l + o(l)), 
quand t — * 0. 

Notons (ic(i, 0, 0) ^ a courbe hamiltonienne issue de (0,0 avec p(0, = ^o- La 
définition du champ hamiltonien donne x(t, = + i 2 5(l) et, d'après la discussion 
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précédente, (t, £) — > x(i,Ç) est un difféomorphisme pour t assez petit mais disjoint de 0. De 
plus, son jacobien vérifie : 

= ft n -\i + o{t)). 



(3.50) 



dx 



dit,® 



3.5. Calcul pour les temps intermédiaires. 

Dans cette partie, on calcule 

(3.51) ^ h = lIo e ~ U{P ~ E)/h K P °ï Sh d *> 

micro-localement près de po G p~ 1 (Eo). Autrement dit, on calcule Op(q)ûh, pour q un symbole 
à support compact proche de po- 



Soient t\, . . . , ij<- les temps strictement positifs définies dans (12.5p . On se donne une fonction 
X G Cq°(M.) égale à 1 près de et à support assez proche de (et donc x es t nulle près des 
tk). Soit x S Cq°(]0, +oo[) une fonction comme dans la partie précédente, à support assez 
proche de et vérifiant x(t — u) = 1 près de pour un certain v > petit. Il existe alors une 
fonction x £ Co°(lR), nulle près de et des tk, telle que, pour t G [0, T], 

K 

(3.52) + "-**) + X(*) = l- 

fe=l 

Ainsi, ïï^ s'écrit 

(3.53) û h = B 1 + B 2 + ~ F x(t)e- U{P - E)/h f(Po)S h dt, 

h Jo 

avec 

x . 

(3.54) fî 2 := £ 1 / x(t - v - t k )e- it ( p - E y h f(P )S h dt. 

k=i •* 



Rappelons que, d'après [21], l'évolution est un opérateur intégral de Fourier d'ordre : 
(3.55) e-^-^G^At), 
dont la relation canonique est donnée par le flot hamiltonien : 

rn , A t = {{x, £, y, v); (x, = K t (y, v)} 

(3.5b) , TT ., . 

La notation cl signifie que le symbole de l'OIF à un développement en puissance de h. De 
plus, le symbole et chaque terme de son développement sont holomorphes par rapport à F\. 

Comme le micro-support de f(Po)Sh est contenu dans {(0, £); Ç 2 /2 G supp/}, le théorème 
d'Egorov implique que le dernier terme dans (|3.53|) est micro-localement nul près de pq. En 
effet, pour t G supp 1[o,t]X(*)> l'opérateur dévolution e - lt ( p - E )/ h envoie le micro-support de 
f(Po)Sfi hors du support de q. Rappelons que "A = B micro-localement près de po" signifie 
que Op(q)A = Op(q)B + 0(h°°). Comme B\ est bien décrit par les sections 1331 et 13.41 il nous 
reste à étudier B2. 
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On peut écrire 

(3.57) B 2 = J2 e -^ k )(P-E)/h^ [ + °° m e- lt ( p - E)/h f(P )S h dt). 
D'après la partie précédente, on sait que 

(3.58) - / x(t)e- u{p - E)/h f{P )S h dt = h~ïb(x, h)é^ x)/h + 0(h°°), 

où b(x, h) est un symbole classique à support dans un compact proche de 0. En particulier, 

(3.59) b(x, h)é^' h G v° cl (À, b ,n), 

où Q C M. 2n est un petit voisinage de {(0,£); Ç 2 = 2(E - V(0))}. La notation V%{A,a ,U) 
désigne l'ensemble des distributions lagrangiennes semi-classiques, u, d'ordre a, de variété 
lagrangienne A, de symbole principal ao (défini sur le fibré de Maslov), de symbole classique 
(c'est à dire que dans toute représentation à l'aide d'une phase non-dégénérée, le symbole de 
l'intégrale oscillante a(x,9,h) vérifie a ~ Ylj>o a j( x >@)h J ) et de micro-support inclus dans 
U C M 2n (c'est à dire Op(l - q)u = 0(h°°) pour tout q G S(l) avec q = 1 près de U). Pour 
toute les questions concernant les distributions lagrangiennes, on renvoie à la section 1.2.1 de 
|18j pour un exposé semi-classique et à |15j dans le cas classique. 

Mais alors, d'après (|3.55p . (|3.56j) . (|3.59j) . le théorème d'Egorov et le théorème 1.2.12 (iii) 
de |18j sur le calcul de l'application d'un OIF à une intégrale oscillante, on a 

e -i(u+t h )(P-E)/h h -$ b fa h)e m*)/h G h -kv° cl (exp({iy + t k )R p )Â, b k ,exp((v + t k )B p )(Q)) . 
En notant A k := exp((f + t k )H p )A C A et £l k := exp((i^ + t k )ïlp)(Çl), on a donc, 

K 

(3.60) B 2 = Y,B% avec B k 2 G hr\v\[k k , b k , tl k ), 

k=0 

micro- localement près de p$. Les A k coïncident avec ceux définis en (|2.6p et les symboles des 
distributions lagrangiennes ainsi que les 6q sont réguliers par rapport à F\. 

Finalement remarquons que, micro- localement près de poi on a 

(P - E)B k =(P -E)'-J X{t~y- t k )e-^ p ~ E y h f{P )S h dt 
= - ! x(t-^-t k )dte-^ p - E ^ h f(P )S h dt 

(3.61) = f d t (x)(t -v- t k )e-^ p - E ^ h f(P )S h dt. 

JR 

En utilisant le théorème d'Egorov et les propriétés du support de dt(x), il vient 

(3.62) (P-E)B 2 k = 0, 
micro- localement près de pç. 
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3.6. Convergence vers une mesure de Wigner. 

Proposition 3.6. — II existe une mesure de Radon positive fi h -i sur M 2n telle que, pour 
tout q G C£°(R 2n ), 

(3.63) h(Op(q)u h ,u h ) — ► (q,(i h -i), 

quand h — > 0. 

Démonstration. — Par linéarité, on peut supposer que q est à support compact près d'un 
point pq. Soit v(x) G C^°(R n ) avec v = 1 près de et vg(x) = v(x/ô) pour 5 > 0. D'après 



l'hypothèse (H4)' et comme l'intersection de deux variétés Aj et est forcément un fermé, 
la régularité de la mesure combinée avec le lemme d'Urysohn implique qu'il existe une suite 
de fonctions f m {x,^) G C^°(M 2n , [0, 1]) telles que, sur le support de q, 

(3.64) mes„(supp(l - f m ) D A) — ► 0, 
lorsque m — ► +oo et A restreint au support de f m ne se recoupe pas. 

Plus explicitement, les f m sont construites de la façon suivante : comme le fermé Aj D A^ 
est de mesure nulle dans la variété lisse Aj, la régularité de la mesure entraine l'existence 
d'un ouvert UJ 1 de Aj, de mesure inférieure à 1/m, qui contient Aj n A&. Comme UJ 1 est 
la restriction à Aj d'un ouvert V™ de M 2 ™, le lemme d'Urysohn nous fournit une fonction 
fj 1 G Cg°(]R 2n ), de support inclus dans VJ 71 , égale à 1 au voisinage de Aj n A&. De la même 
manière, on construit f™ en intervertissant les indices j et k. On pose alors f m = 1 — /j"/™ 
que l'on tronque hors du support de q. 

Soit q G C^°(M 2n ) à support proche de celui de q et égal à 1 près de celui-ci. D'après (|3.9p . 
(|3.53[) . la discussion après (|3.56|) et les propriétés du calcul pseudo-différentiel, on a 

(3.65) Op(q)u h = Op(g)(5i + B 2 ) + R et Op($H = Op{q){B 1 + B 2 ) + Ro, 

avec R ,Ro = 0~^(ï) pour tout [3 > 0. Le lemme [331 (13^71) et (pTëUD et la proposition 1.2.2 
de [18 1 sur la norme L 2 d 'une intégrale oscillante impliquent 

(3.66) v&Bi = Ri = 0(Vôh~^) 

(3.67) Op(l - / m )(l - v s )B 1 =R 2 = om^ih-ï) + O m (hh) 

Op(f m )(l-v s )B 1 = h-^BT' S 

(3.68) avec B™' S G V° cl (À, (1 - v s )f m \ Ah bo, fi n supp/ m (l - v s )) 

(3.69) Op(l - f m )B 2 = R 3 = om^ih-h) + O m (hh) 



Op(f m )B 2 = h-ïJ2B, 



K 

2 

k=0 



(3.70) avec B*> m G V° cl (A k , b k , ft k D supp f m ) 

Ici, une fonction est un cv^—x^^l) si elle tend vers lorsque m — ► 00. De même, une fonction 
est un O m {h) si, pour tout m fixé, elle est un 0(h). En particulier, 

(3.71) Op(?H = 0-P{h-*), 
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pour tous p > et q G C^(M. 2n ). 
On a alors 

(Op(q)u h ,u h ) =(Op(<zK,Op(gH) +0(h°°) 

= (R + Op{q){Ri + R 2 + R 3 ), Op(g)« ft ) 

+ (Op(g)(i?r' 5 + E fi 2' m )^o + Op(ç)(i2i + iî 2 + iîs)) 

+ / i - 1 (0 P ( (? )«' 5 + E B 2' m ), Op®^ + E B 2' m )> + ^°°) 

=/ l - i (o P (g)(i?r' 5 + e B k 2 n, (K' 5 + E B 2 m )) 

(3.72) + /i~ 1 (o (l) + 0{V6) + O s , m (h*)). 

Comme les K + 1 variétés lagrangiennes A et A& ne s'intersectent pas sur fî n supp / m (l — vg) 
et fifc H supp/ m , les distributions lagrangiennes B™' S et B 2 ' m ont des micro-supports deux à 
deux disjoints. Ceci entraine 

(3.73) (Op(q)BT' S , B k ' m ) = O m 4h°°) 

(3.74) (Op(q)B k ' m , B k ' m ) = O m {h°°) pour k^k. 

Le calcul de l'application d'un opérateur pseudo-différentiel sur une distribution lagran- 
gienne (théorème 1.2.8 de [18]) donne 

(3.75) Op(g)Dâ(A, a , U) C X>â(A, g, A x a , U). 
En particulier, 

o v { q )B^ 5 g 2?â(Â )g (i-^)/ m | Aj| 6o,nnsupp/ m (i-«i)) 

0p(ç)5 2 fe ' m G 2?â(A fc , ç/ m | A)fe 6 Q> n fc n su PP / m ). 
Un calcul évident montre alors que 

(3.76) (Op(q)B™' 5 ,B™' 5 )= fq\(l-v s )f m b \ 2 d l X + O m>s (h) 

(3.77) (Op(q)B k 2 > m ,B k > m ) = [ q\f m b k \ 2 df, + O m (h), 

où // est la mesure de Lebesgue sur À ou Aj. Comme 6q est régulière, on a 

(3.78) (Op(q)B k > m ,B k > m }= [ ç|6g| 2 ^ + o m ^ 00 (l) + O m (/ i ). 

JAk 

Pour la partie concernant B™' S , le symbole bo donné par (|3.48p n'est pas régulier en car 
| det <p%t 1 2 explose quand x — > 0. Mais en effectuant le changement de variable dont le Jacobien 
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est donné en (|3.50p et en utilisant l'estimation (|3.49p . on peut écrire 

dx 



(Op(q)B?' S ,B?' 5 )=J g|(l -«,)/, 



m| 2 |&o| 2 



d(t,0 



(3.79) 



q\b \ 2 (x(t,0y^(x(t,0)) 



dtdÇ + O m>s (h) 
dx 



d(t,0 

l) + W (l)+Om#)- 



dtdÇ 



Dans la formule précédente, l'intégration se fait sur [0, £o[ x y^C^ ~~ y(0))S™~ 1 pour Eq assez 
petit. On reconnaît une mesure dans le premier terme de (13.79p . La proposition résulte alors 
de (|3.72|) . (|3.73p . (|3.74p . (|3.78|) et (|3.79|) : on fixe d'abord ô assez petit et m assez grand et 
finalement, on fait tendre h vers 0. 

Ceci conclut la preuve pour E = Eq + hF\. Dans le cas E = Eq + hE\ + Oh—,o(h), il résulte 
des constructions précédentes que les symboles b et b$ qui apparaissent dans (|3.78|) et (|3.79|) 
sont holomorphes par rapport a F\ et que les termes d'erreur sont localement uniformes par 
rapport â F\. Il suffit alors de prendre F% = E\ + o^,_»o(l) et d'utiliser la formule de Taylor 
eni^.Ei. □ 



3.7. Calcul de la mesure fi^-i. 

Dans cette partie on calcule la mesure de Wigner p h -i définit dans la proposition 13.61 On 
a la proposition suivante 

Proposition 3.7. — La mesure p h -i possède les propriétés suivantes qui la caractérisent 
uniquement : 

i) La mesure p, h -i est supportée dans la surface d'énergie Eq : 

(3.80) supp^-iCp- 1 ^) 
H) La mesure fi^-i est nulle dans la zone entrante : 

(3.81) p h -i = dans {(x,£); \x\ grand et cos(x,Ç) < —1/2}. 
iii) La mesure fi^-i vérifie l'équation d'évolution : 

(R p + 2ïmEi)iJ, h -i =(^-Vy(i)9 c + 2ImEi)/i H 



(3.82) 



(27ry- n \S(trô x=0 ô le+v{0> 



E 



Démonstration. — L'estimation (13. 2p combinée avec l'existence de la mesure de Wigner (pro- 
position [3TB]) implique (|3.80p . 



La preuve de (I3.8ip résulte de l'existence de la mesure et du fait que Uh = O a (h°°) dans 
la zone entrante (équation (|3.6p et proposition 13.41) . 



Soit q G C^°(M. 2n ) localisé près de po = (xq,Ço)- Comme (P — E)uh = Sh, on a (P 



E)iih = 0(h°°) micro-localement près de po, si xq ^ 0. D'après (H6) E — E = 2iIm.E\h + 
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o(h). En remarquant que H p (g) est le symbole principal de l'opérateur pseudo-différentiel 
i/i~ 1 [P,Op(g)], il vient 

{q, (H p + 21mEi)n h -i)=-{U p (q)-2qTmE 1 ,ii h -i) 
(3.83) = - lim <i((P - P)* Op(<z) - Op(g)(P - £7))u h ,u A ) + o(l) = 0. 

Ainsi, (H p + 2ImPi) ^-î = pour x ^ 0. 

On suppose maintenant que xo = 0. Près de pq, la mesure Ph- 1 est donnée comme somme 
des termes dominants de flgTTgD et (pTf9l) . D'après ([3T62]) . B\ vérifie (P - P)Pf = micro- 
localement près de po et la preuve de (|3.83p montre que ces termes ne contribuent pas dans 
le calcul de (Hp + 2ImPi) p^-i. Donc, 



(q, (H p + 21m E{) p h -x) = - J (H p (ç) - 2q Im E x ) |6 | 2 (^(i, £), V^(z(t, 0)) 
(3.84) = - / (%(p(t,e))-2g(p(t,e))Im£; 1 )|6o| 2 (x(t,e)) 



d(t,o 

dx 



dtd£ 
dtdÇ, 



où p{t,£,) est la courbe hamiltonienne avec donnée initiale (0,£) en t = 0. Comme avant, 
l'intégration se fait sur [0, e [x y / 2(E - V(0))S n - 1 . D'après (|Q7)> . (pOHL (EQ9Î) et f[375Ô]) . 
on a 



(3.85) 



|èo| 2 (x(t,Ç)) 



(27r) i -»|flf(Or(l + 0|^o(l)). 



Soit x S Cq°(M, [0, 1]) avec x — 1 P res de et on note x m (i) = x(* m )- On voit Xm comme 
une fonction définie sur un voisinage de po S M? n en posant Xm(x,v) = XmW si x = x(t,£) 
(rappelons que (t, £) — > x(t, £) est le changement de variables défini avant (|3.50p ). En utilisant 
la discussion précédente, il vient 

{q,(H p + 2ImE 1 )p: h -x) 

=(x m q, (H p + 2lmE 1 ) p h -i) + ((1 - Xm)q, ('H v + 21mE 1 )fi h -i) 
= - ((Hp - 2ImP 1 ) ( Xm q),p h -i) + 

{x(tm)(d t q-2qlmE 1 )+mxXtm)q)(2Tr) 1 - n \S(0\ 2 (l+o t ^ (l))dtdÇ 



(3.86) 



Oim- 1 ) - JJ^ mx'(tm)q(p(0, 0) (2vr) x — | 1 2 (1 + o^ (l)) dt d£ 
g(0,0(27r) 1 -"|%)| 2 ^ + 0^00(1) 
g(0,0(27r) 1 -"|%)| 2 ^, 



où l'intégration se fait sur y / 2(E - V r (0))S n_1 . 



□ 
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4. Contre-exemple quand l'hypothèse |(H4)| n'est pas vérifiée. 



On donne ici un exemple (en dimension n = 1 puis en dimension quelconque) de potentiel 
et de fonction S qui ne vérifient pas |(H4)"| et dont la mesure limite n'est pas unique. Pour fixer 
les idées, on suppose que E = Eq. 

En dimension n = 1, considérons un potentiel satisfaisant les hypothèses |(H1)[ |(H2)| et 
|(H3)| de la forme suivante : 

A V(x) 








Figure 3. Le potentiel V(x). 

L'ensemble A ressemble alors à : 

A t 




Figure 4. La variété A avec A_ = A + dans la zone quadrillée. 



Soit S une fonction vérifiant |(H5)| D'après la partie [331 la solution u^, micro-localement 
dans la zone hachurée, est de la forme 

(4.1) u h = a-{x,h)é^- {x)/h + 0(h°°), 

où a~ est un symbole classique et ip = £ dx = ^(s) 2 ds. En particulier, si le support 

de g G Cg°(M 2n ) est contenu dans la zone hachurée, 

(4.2) h(Op(q)u h: u h ) — ► (q,fj, h -i), pour h -v 0, 
et fih-i vérifie les conclusions du théorème 12.11 
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Par contre, micro-localement dans la zone quadrillée, Uh est de la forme 

(4.3) u h = a-{x,h)e ilf '~ (x)/h + a + (x,h)é i ' +(x)/h + 0(h°°), 

où a + et vérifient le même type de propriétés que a~ et ip- . En particulier, 

(4.4) ip-{x) = \ Çdx= Çdx + / Çdx = A + ip+(x), 

Jj-(x) y 7 _(o) Ji+(x) 

où A = Ç(s) 2 ds > est l'action associée à la courbe 7-(0). En outre, en utilisant (|3.48|) . 
le fait que les a* vérifient la même équation de transport dans (la projection en x de) la zone 
quadrillée et la théorie de Maslov, on trouve 



où v est l'indice de Maslov associé à la courbe 7-(0) et r, r G]0, +oo[, € M sont des constantes. 
Pour que Oq et soient non nuls, on suppose que les transformées de Fourier de S qui 
apparaissent dans (|4.5p sont non nulles. Ainsi, si le support de q G C^°(lR 2n ) est contenu dans 
la zone quadrillée, 



h(Op(q)u h ,u h ) = j g(x,VV.(x))(|aô| 2 + |a+| 2 + 2Re (aô a+ e iA ' h )) dx + o h ^ (l) 

(4.6) = J g(x,VV'.(x))(|aô| 2 + |a ( [| 2 + 2|a "||a ( [|cos(6l + ^l//i)) ( ia; + o/ l ^o(l)- 
Comme A ^ 0, on obtient facilement la 

Proposition 4-1- — Sous les hypothèses précédentes et pour tout v G [—1,1], on peut 
trouver un ensemble de h adhérant à tel que 

(4.7) h{Op(q)u h ,u h ) — ► y q(x, VV>.(x))(|aô| 2 + |a^| 2 + 2^|oq | |) dx, 
quand h tend vers dans cet ensemble. 

En particulier, la mesure limite n'est plus unique. La mesure du théorème 12.11 correspond 
à v = 0. Il se peut que la mesure soient nulle en dehors d'un compact. En dimension plus 
grande, en s'inspirant de [3], on peut construire un potentiel tel que localement près d'un 
point po l'ensemble A soit la réunion de 2 variétés A_ et A + qui coïncident (figure [5]). 

Comme en dimension n = 1, micro-localement près de po, la fonction Uh s'écrit 

(4.8) u h = a-(x,h)e i ' p ~ ix)/h + a + (x,h)e li,+ ix)/h + 0(h°°), 

La différence des phases A(x) = ip-(x)—ip + (x) est encore l'action le long de la courbe 7_ i3; (0) : 



(4.9) A(x) = Ç-dx. 

Mais on remarque que cette action ne dépend pas de x. En effet, si 7_ jXl (0) et 7_ jX2 (0) sont 
deux chemins, le fait que les l, soient dans {x = 0}, la formule de Stokes et le caractère 
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x = 



Figure 5. Un exemple de potentiel, égal à l'oscillateur harmonique dans un petit 
secteur angulaire, qui re-focalise en tous les rayons issus d'un secteur angulaire en 
0. L'ensemble A, avec A_ = A + près de po, est décrit à droite. 



lagrangien de A impliquent 



£ • dx 



(4.10) 



Ç • dx 



7-,«a(0) 



/ £ • dx + / 

J'y-.*, (0) Jlo 



i-dx+ I Ç-dx- l 
h-,xx(o) Jl 2 J-y 



7-,«a(°) 



£ ■ dx + I £ • dx 

h 



dÇ A dx = 0. 



On peut aussi dire que comme A_ et A + coïncident, on a V^_ = et donc ip- et V+ 

diffèrent d'une constante. Cette remarque permet d'adapter la preuve de la proposition 14.11 à 
cette situation. 



5. Preuve de la proposition 12.31 



En notant a#6 le symbole de l'opérateur Op a o Op b, on a 
q(p - E ) 2 =(p - £ ) (?#(P " Eq)) + \.{V ~ £o)H p q + ^ 2 r(x, £, h) 

= (p - E )#(q#(p - E )) - ^R p (q#(p - Eq)) + ^{p - E )R P q + h 2 r(x, Ç, h) 



(5.1) 



= (p - E )#{q#(p - Eq)) + h 2 r(x, Ç, h), 



où r G 5 c i(l) change de lignes en lignes et à son support inclus dans le support (compact) de 
q, modulo S(h 00 (x)^ 00 (Ç i )^ 00 ). Mais alors, en utilisant que Uh = 0(h~ 1 / 2 ) micro-localement 
près du support de q, \\Sh\\ = 0(1) et (Op(q)uf l ,Uh) = 0{h~ l ) (d'après le théorème 12. ip . il 
vient 

<Op(g(p - E ) 2 )u h ,u h ) =(Op(p - Eq) Op(q) Op(p - E )u h , u h ) + O(h) 

=(Op(<?) Op(p - E)u h , Op(p - E)u h ) + OpOzHIIHShll + 0(h) 



(5.2) 



=<Op(g)5 h) 5 fc ) + 0(^ 
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Grâce à la formule sur l'application d'un opérateur pseudo-différentiel sur une distribution 
lagrangienne (théorème 1.2.8 de |18j) ou en faisant le calcul explicitement, on obtient 

Op(q)S h =Op(q)g(x)S h + 0(h 00 ) 



= (2vr)- n /i- n / 2 j O V (q)g(x)e lx < /h S(0 + 0(h°°) 

(5.3) ={2-K)- n h- n ' 2 j q(x, 09(x)e iX</h S(0 dÇ + 0(h), 
pour g{x) G C^iW 1 ) avec g = 1 près de 0. En particulier, (|5.2p devient 
(Op(q(p-E ) 2 )u h ,u h ) = (2<K)- 2n h- n JJ j q(x,0g(x)e ix </ h S(0e- ix - e / h ~S(e)dxdÇd9 + O(h). 
La méthode de la phase stationnaire en variables (x, 9) donne alors 

(5.4) (Op(ç(p - E ) 2 )u h , u h ) = (2tt)-" J q(Q, 0\S(i)\ 2 + D(h), 
ce qui démontre la proposition. 



6. Equivalence entre les hypothèses |(H4)| et (H4)'[ , 



Lemme 6.1. — Soient Ai et A2 deux sous-variétés de M n de même dimension m. Alors 
l'ensemble fermé Ai n A2 est de mesure nulle dans Ai si et seulement si il est de mesure nulle 
dans A2. 

Démonstration. — On travaille localement près de po G Ai n A2. Il est possible de trouver, 
E, un sous-espace vectoriel de W 1 de dimension m — n qui soit supplémentaire à T po Ai et 
à T po A2. Pour cela, on peut utiliser la propriété suivante : si A, B sont deux sous-espaces 
vectoriels d'un espace vectoriel G, alors A U B = G si est seulement si A = G ou B = G. 

On note F = Tp Ai et II la projection linéaire de W 1 sur F parallèlement à E. Alors, II est 
un difféomorphisme local de Ai sur F et de A2 sur F. De plus, les mesures de Lebesgue sur A, 
sont transportées en des mesures sur F. Ces deux mesures ainsi que la mesure de Lebesgue 
sur F sont absolument continues les unes par rapport aux autres. En particulier, Ai n A2 est 
de mesure nulle dans Ai ssi II(Ai n A2) est de mesure nulle dans F pour une des trois mesures 
précédentes ssi Ai n A2 est de mesure nulle dans A2. □ 

On rappelle que si p(t) = (x(t),Ç(t)) est une courbe hamiltonienne dans p~ 1 (Eq), alors 
p(t) = (x(—t), —£(—£)) en est aussi une. Donc, 

G V2(£o - nO))^" 1 ; 3t < IL, exp(tH p )(0, £) = 0} 

= -{^ G y/2{Eo - y(0))S' 1 - 1 ; 3t > U x exp(tH p )(0, £) = 0}. 



En particulier, dans l'hypothèse (H4) , on peut remplacer t > par t 7^ 



Lemme 6.2. — Supposons (H1)--(H3). Les hypothèses |(H4)| e£ |(H4)' sont équivalentes. 



L'équivalence entre les hypothèses |(H7) | et | (H7) ' 



se démontre de la même façon. 
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Démonstration. — Supposons | (H4) | et considérons l'intersection de Aj et avec j < k près 
d'un point po. La variété Aj est paramétrée par l'application p(t,£) = exp(iH p )(0, £) avec t 
au voisinage de tj et £ G ^/2{Eq — V(0))§ n_1 au voisinage de Çj. 

Donc, un point p = p(t, £) est dans Aj n ssi il existe s proche de t k tel que p(s, £) = p(t,Ç) 
ssi il existe u = s — t proche de tu — tj tel que p(u, £) = (0, £) ssi 



£ G M i)È = G ^2(^0 - ^(0))S™- 1 proche de fa 

3t>0 proche de tk - tj H x exp(tH p )(0, £) = 0}. 
En utilisant le théorème de Fubini pour séparer les variables t et £, on a 

mes n Aj H A k = jT l AjnAfc c^ < ^ l p (i,ç) eAj l^ eMj:k dt 
< y IfeMf, k d£ = mes n _i M jtk . 



(6.1) 

Comme mes^—iMj^ 



d'après (H4) l'hypothèse (H4)' est vérifiée. 



Réciproquement, supposons (H4)' Soit Ço £ \/2(-Eb — V(0))S n . En utilisant le même 
raisonnement que (|6.ip mais en sens inverse, on trouve 



mes n _i{Ç G ^2(^0 - ^(O))^"" 1 proche de Ç ; 3t > n x exp(tH p )(0, £) = 0} 
= ^ M , fc < ^ mes„, Â n Afc. 

D'après l'hypothèse (H4)', cette quantité est nulle. L'hypothèse |(H4)"] découle alors d'un ar- 
gument de compacité en £o- D 
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